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ABSTRACT 

We identify the largest known lensed images of a single spiral galaxy, lying close to the centre of the 
distant cluster MACS J1149. 5+2223 (z = 0.544). These images cover a total area of ~ 150 □" and 
are magnified 200 times. Unusually, there is very little image distortion implying the central mass 
distribution is almost uniform over a wide area (r ~ 200 kpc) with a surface density equal to the critical 
density for lensing, corresponding to maximal lens magnification. Many fainter multiply-lensed galaxies 
are also uncovered by our model, outlining a very large tangential critical curve, of radius r ~ 170 kpc, 
posing a potential challenge for the standard LCDM-Cosmology. Because of the uniform central mass 
distribution a particularly clean measurement of the mass of the brightest cluster galaxy is possible 
here, for which we infer stars contribute most of the mass within a limiting radius of 30 kpc, with a 
mass-to-light ratio of M/Lb — 4:.5{M/ L)q. This cluster with its uniform and central mass distribution 
acts analogously to a regular magnifying glass, converging light without distorting the images, resulting 
in the most powerful lens yet discovered for accessing the faint high-z Universe. 

Subject headings: gravitational lensing , galaxies: clusters: individual: MACS J1149. 5-1-2223 , dark 
matter 



1. INTRODUCTION 

Multiply-lensed images of distant galaxies are commonly 
seen near the centres of distant galaxy clusters. Thin tan- 
gentially stretched images, including giant arcs, follow an 
approximate "Einstein ring" , interior to which radially di- 
rected arcs are sometimes found, pointing towards the cen- 
tre of mass. The formation of such strongly lensed images 
requires the projected mass density to exceed a critical 
value, Scrit, inside the Einstein radius, given by funda- 
mental constants and the angular-diameter distance of the 
lens, di, and of the source d^, and their separation, dis, 
such that Scrit = i^IvD^- This mass density is a little 
less than 1 g/crn^ for lenses at intermediate distances and 
comprises mainly dark matter (DM) whose nature is still 
unknown, with only a ~ 5-15% contribution from gaseous 
baryonic material (Biviano & Salucci, 2006, Umetsu et al., 
2009). 

Fine examples of lensing by galaxy clusters are found at 
intermediate redshifts, including the most distant cluster 
discovered by Zwicky, C10024-I-17 {z = 0.39), and one of 
the richest clusters discovered by Abell, A1689 {z — 0.18), 
with the largest known Einstein ring, r ~ 45" (Broad- 
hurst et al., 2005). For such clusters many tens of mul- 
tiply lensed images have been identified in deep Hubble 
images (Broadhurst et al., 2005, Halkola ct al., 2008, 
Limousin et al., 2008, Zitrin et al., 2009a, b), leading to 
accurately measured central surface mass distributions. 
The Einstein radii of these massive clusters are found to 
be larger than predicted in the context of the standard 
LCDM cosmological model (Broadhurst & Barkana, 2008, 
Sadeh & Rephaeli, 2008), based on the "Millenium" sim- 
ulation (Springel et al., 2005). This discrepancy is em- 
pirically supported by the surprisingly concentrated mass 
profiles measured for such clusters, when combining the 
inner strong lensing with the outer weak lensing signal 
(Gavazzi et al., 2003, Broadhurst et al., 2008, Limousin ct 



al., 2008, Umetsu et al., 2009, Zitrin et al., 2009a, Oguri 
et al., 2009, Donnarumma et al., 2009) boosting the crit- 
ical radius at fixed virial mass. The abundance of giant 
arcs may help constrain the total lensing cross section and 
is variously claimed to be at odds with standard LCDM 
(Bartelmann et al., 1998), though recent work favours the 
consensus (Dalai et al., 2004, Wambsganss et al., 2004, 
Horesh et al., 2005, Hennawi et al., 2007) here the selec- 
tion effects are considerable and more thorough surveys 
should help (Hennawi et al. 2007). 

The magnification generated by massive clusters has 
consistently led to the discovery of the highest redshift 
galaxies (Ebbels et al., 1995, Franx et al., 1997, Frye & 
Broadhurst, 1998, Bouwens et al., 2004, Kneib et al., 2004, 
Zheng et al., 2009), with the current record standing at 
z ~ 7.6 for a galaxy behind A1689 (Bradley et al., 2008) 
and magnified by nearly a factor of ~ 10. Although lens 
magnification, /i, reduces the accessible area of the source 
plane by l//x, it enhances the flux of faint galaxies by 
fi, with a net positive effect for the most distant galax- 
ies lying on the steep exponential tail of the luminosity 
function (Broadhurst, Taylor & Peacock, 1995, Bradley et 
al., 2008). Lensing provides additional spatial resolution 
by stretching images, producing spatially resolved details 
and evidence that outflowing gas is widespread at z > 4 
(Frye, Broadhurst & Benftez, 2002). 

With the goal of discovering high redshift galaxies and 
to better define the mass profiles of galaxy clusters in gen- 
eral, we have combined data for a sample of well studied 
clusters. In this process we have uncovered the unusual 
lensing properties of MACS J1149. 5-1-2223 (z=0.544), a 
cluster originally identified in the highly complete sam- 
ple of the most X-ray luminous clusters in the Universe 
(Ebeling et al., 2007). In §2 we describe the observations, 
in §3 we describe the lensing analysis and results, in §4 
we address the BCG, and in §5 discuss our conclusions. 
Throughout the paper we adopt the standard cosmology 
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(ilnio = 0.3, flAo = 0.7, h = 0.7). Accordingly, one arc- 
second corresponds to 6.4 kpc/hjo at the redshift of this 
cluster. The reference center of our analysis is fixed at 
the center of the cD galaxy: RA = 11:49:35.70, Dec = 
+22:23:54.8 (J2000.0). 

2. OBSERVATIONS 

The central region of the luminous X-ray cluster MACS 
Jl 149. 5-f 2223 (z = 0.544) has been imaged in April 2004 
and in May 2006, with the Wide Field Channel (WFC) of 
the ACS installed on HST. Integration times of ^ 4500s 
were obtained through the F555W and the F814W filters. 
We retrieved these images (PI: Ebeling, proposal ID: 9722) 
found in the Hubble Legacy Archive. Several large blue 
spiral galaxy images are clearly visible near the central 
brightest cluster galaxy (Figures 1-3). On closer inspec- 
tion, individual HII regions and spiral arms are repeated 
very clearly in images 1.1 and 1.2, demonstrating beyond 
question that these are images of the same source, even 
though they do not appear as thin distorted arcs. The 
other central spiral galaxy images seen in Figure 1-3 are 
also images of the same source but with differing mirror 
symmetry (parity), as we show below. Image 1.2 is the 
largest, covering an area of ~ 55 □", and in total > 150 □" 
is subtended by all the images of this source, several times 
greater than the largest giant arc known (Soucail et al., 
1987). 

Many other faint lensed galaxies are also visible, gener- 
ally at larger distances from the centre (marked in Figure 
2) indicating that they lie at higher redshift and most of 
which we have been able to securely identify as sets of 
multiply-lensed background galaxies as detailed below. 
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Fig. 1.— The central region of the cluster MACS J1149.5-I-2223 
(z = 0.544) imaged with Hubble/ACS in V and I-bands, showing 
the several large, well resolved, blue spiral galaxy images centered 
on the brightest red cluster member galaxy (cD). The small de- 
lensed source, zoomed-in x4 to show internal details, is inset near 
its predicted location demonstrating the very large magnifications 
generated by the central gravitational field of this cluster. 



3. LENSING ANALYSIS 

We apply our well tested approach to lens modelling, 
which we have applied successfully to A1689 and C10024, 
uncovering unprecedentedly large numbers of multiply- 
lensed images (Broadhurst et al., 2005, Zitrin et al., 
2009a). The full details of this approach can be found 
in these papers. Briefly, the basic assumption adopted is 
that mass approximately traces light, so that the photome- 
try of the red cluster member galaxies is the starting point 
for our model. 

Cluster member galaxies are identified as lying close to 
the cluster sequence by the photometry provided in the 
Hubble Legacy Archive. We approximate the large scale 
distribution of matter by assigning a power-law mass pro- 
file to each galaxy, the sum of which is then smoothed. 
The degree of smoothing and the index of the power-law 
are the most important free parameters. A worthwhile 
improvement in fitting the location of the lensed images 
is generally found by expanding to first order the gravi- 
tational potential of the smooth component, equivalent to 
a coherent shear, where the direction of the shear and its 
amplitude are free, allowing for some fiexibility in the rela- 
tion between the distribution of DM and the distribution 
of galaxies which cannot be expected to trace each other 
in detail. The total deflection field ariO), consists of the 
galaxy component, dgai{0), scaled by a factor Kgai, the 
cluster DM component doMid), scaled by {1-Kgai), and 
the external shear component dex {&) ■ 

driO) = Kgaldgalie) + (1 - Kgal)dDM{0) + dex{0), (1) 

where the deflection field at position 9m due to the exter- 
nal shear, dex{Om) = {aex,x,aex,y), is given by: 

aexA^^m) = \l\ cos{2(j)j)Axm + 111 sm{2(j)^)Aym, (2) 
aex,yiOm) = ItI sin(20-y) Ax^ - I7I cos(2(?!)-y)Ay™, (3) 
where (Axm, Aym) is the displacement vector of the 
position 9m with respect to a flducial reference position, 
which we take as the lower-left pixel position (1, 1), and 
(f)j is the position angle of the spin-2 external gravitational 
shear measured anti-clockwise from the x-axis. The nor- 
malisation of the model and the relative scaling of the 
smooth DM component versus the galaxy contribution 
brings the total number of free parameters in the model 
to 6. 

We lens all well detected candidate lensed galaxies back 
to the source plane using the derived deflection field, and 
then rclcns this source plane to predict the detailed ap- 
pearance and location of additional counter images, which 
may then be identified in the data by morphology, inter- 
nal structure and colour. The fit is assessed by the RMS 
uncertainty in the image plane: 

RMSfmages = " ^if + {V^ " y^f) I N.mages , (4) 

i 

where x^ and are the locations given by the model, and 
Xi and yi are the real images location, and the sum is 
over all Nimages imagcs. The best-fit solution is unique in 
this context, and the uncertainties are determined by the 
location of predicted images in the image plane. 

Importantly, this image-plane minimisation does not 
suffer from the well known bias involved with source plane 
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minimization, where solutions are biassed by minimal scat- 
ter towards shallow mass profiles with correspondingly 
higher magnification. The model is successively refined as 
additional sets of multiple images are identified and then 
incorporated to improve the model (Zitrin et al., 2009a). 




Fig. 2. — Large scale view of the multiply lensed galaxies identi- 
fied by our model. In addition to the large spiral galaxy system 1, 
many other fainter sets of multiply lensed galaxies are uncovered by 
our model. The white curve overlaid shows the tangential critical 
curve corresponding to the distance of system 1. The larger crit- 
ical curve overlaid in blue corresponds to the average distance of 
the fainter systems, passing through close pairs of lensed images in 
systems 2 and 3. This large scale elongated "Einstein ring" encloses 
a very large critically lensed region equivalent to 170 kpc in radius. 
For this cluster one arcsecond corresponds to 6.4 kpc//i7o, with the 
standard cosmology. 
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Fig. 3. — The observed images are shown with (top row) and 
without (bottom row) the light of the cluster member galaxies, and 
compared with our model generated images (middle row). For the 
model we take as input the pixels belonging to 1.1, the least dis- 
torted and cleanest image in system 1, and delens these pixels back 
to the source plane (see inset in figure 1) and then relens the source 
plane to generate counter images. It is clear that our modelling is 
successful in demonstrating the multiply-lensed origin of all 5 ob- 
served images, corresponding to a single distant spiral galaxy. 
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Fig. 4. — 2D surface mass distribution (k), in units of the critical 
density. Contours are shown in linear units, derived from the mass 
model constrained using 33 multiply lensed images seen in Figure 2. 
Note, the central mass distribution is shallow, and rounder in shape 
than the critical curves. 

The derived surface mass distribution is found to be 
very nearly uniform within the central ~ 200 kpc (Fig- 
ures 4-5), with very little uncertainty, as must be expected 
given the very large and undistorted images observed. The 
value of the uniform surface mass in this central region is 
the critical value for generating multiple images, about ~ 
0.50 g/cm? at our estimated redshift for the source, see 
below. 

The total magnification of the spiral galaxy we derive is 
about, ~ 200, when summed over all five images, forming 
the largest known images of any lensed source, and is inde- 
pendent of the unknown source redshift, and given by the 
ratio of the area of images divided by the area subtended 
by the deprojected source. This can be appreciated in Fig- 
ure 1 where we plot for comparison the unlensed source 
and its modelled location on the sky, whose diameter is 
estimated to be ~ 0.9", typical of unlensed spiral galaxies 
at intermediate redshifts. 

The critical curve corresponding to the spiral system is 
shown in Figure 2, together with the larger Einstein radius 
derived for the fainter galaxies, which lie at higher redshift, 
typically in the range z ^ 1.5 — 2.5, based on other massive 
clusters where redshift measurements have been made or 
photometric redshifts obtained with more colour informa- 
tion (Broadhurst et al., 2005, Zitrin et al., 2009a). From 
this we can roughly calculate the relative lensing distance 
of the spiral system 1, from the ratio dis/ds- For a mean 
background redshift of z = L5, the spiral system is at 
z ~ 1.2, and at the other extreme for a background red- 
shift of 2 — 2.5 the maximum predicted redshift of the 
spiral is z ~ 2. Assuming an average background depth of 



z ~ 2±0.5 yields Zgpirai = l-5_o'3i due to the shape of the 
dis/dg growth as a function of source redshift for a clus- 
ter at z = 0.544. The shape, location and magnification 
of the critical curves are unaffected by this uncertainty. 
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The tangential critical curve is elongated, with a major 
axis of ~ 80" at a redshift of '--^ 2 (Figure 2), reflecting 
the somewhat elongated distribution of matter (Figure 4). 
Note, the mass distribution is more symmetric than the 
shape of the tangential critical curve, the form of which is 
sensitive to asymmetry. The circular equivalent Einstein 
radius contained within the critical area is ~ 27" for faint 
sources at z ~ 2, corresponding to the outer tangential 
critical curve drawn in Figure 2 (light blue curve), which 
in physical units is 170 ±20 kpc at the distance of the clus- 
ter, where the uncertainty is dominated by the uncertain 
mean redshift of the background galaxies. Note that be- 
cause of the weak redshift dependence of lensing distance, 
dis/dg, the uncertain redshifts are only a minor source of 
uncertainty when deriving physical scales. 
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Fig. 5. — Radial surface mass density profile, K{r), in units of the 
critical surface density, i.e. K{r) = T,{r) /T,crit, derived for the range 
of radius covered by all sets of multiple images shown in Figure 2, 
The profile is very flat within 200 kpc, at a level equal to the crit- 
ical density for system 1 (black line). The rise in the centre below 
~ 5" is due to the contribution of the cD galaxy, under which the 
modelled smooth cluster mass component is indicated by the black 
Une. Three geometrically averaged Einstein radii are indicated, cor- 
responding to 16" radius for system 1, and 27" at the mean depth 
of the fainter systems (z ~ 2) and 47" is predicted for the Einstein 
radius at high redshift, 2 ~ 8. 

4. BRIGHTEST CLUSTER GALAXY 

The mass profiles of brightest cluster galaxies are not 
very well constrained presently by lensing, but may help 
elucidate the origin of these poorly understood class of ob- 
jects (von der Linden et al., 2007, and references therein). 
The best studied example with lensing is MS2137 (Gavazzi 
et al., 2003, Donnarumma et al., 2009), where radial 
arcs lie close to this object. For clusters like A1689 and 
C10024-I-17 there are too many massive galaxies in the cen- 
tral region to allow a unique determination of the brightest 
galaxy. Here we are fortunate to have such extended lensed 
images spread over a wide range of radius to help constrain 
the central mass distribution. Furthermore, the contrast 
of the cD mass above the flat cluster profile is clear, al- 
lowing an accurate subtraction of the cluster contribution, 
as shown in Figure 5. We model this galaxy analytically 
as a general power-law mass profile with a projected core 
(Barkana 1998), resulting in a mass of - 1.0±0.2x lO^^M© 
interior to the limiting radius of ~ 30 kpc traced by the 



low surface brightness wings of the light profile (see Figure 
1), where the uncertainly is dominated by the gradient of 
the cD mass profile (Figure 5). 

Integrating the light within this radius yields 2.2 x 
10^^ Lq, in the restframe B-band and hence a mean mass- 
to-light ratio of ~ M/Lb = 4.5±1 {M/L)q. This ratio can 
be fully accounted for by the stars contained in this galaxy, 
for which we obtain M/Lb — 5 {M/L)q, for a single burst 
stellar population formed at z = 3 and viewed at a red- 
shift oi z = 0.544, equivalent to an age of ~ 6 Gyrs, and 
with a mean metallicity of half solar. Other lensing work 
supports relatively low mass-to-light for the cD galaxy in 
MS2137 {M/L < 10 {M/L)q, Gavazzi et al., 2003), with 
less direct estimates relying on dynamical motions show- 
ing evidence for non-stellar DM at larger radius for other 
cD galaxies (Dressier 1979, Gebhardt & Thomas, 2009). 
This demonstrates the need for more thorough studies of 
these enigmatic objects. 

5. DISCUSSION 

The unusually large and undistorted lensed images of 
a spiral galaxy uncovered here requires a nearly uniform 
distribution of matter within the central 200 kpc region 
covered by these images. The formation of multiple im- 
ages requires the value of the central surface density to be 
nearly equal to the critical surface density for lensing. This 
lens corresponds to the case of maximum magnification, 
given by the general relation /X = (((1— (S/Ecrit)^— 7^))~^, 
because the lensing shear, 7, vanishes for a uniform density 
and therefore the magnification diverges when the surface 
density is equal to the critical density, S — "Scrit- These 
conclusions follow from the fundamentals of lensing and 
do not rely on any model. 

In detail we show that the locations and parity of the 
lensed images can be accurately reproduced with a simple 
model where indeed the central mass distribution is nearly 
uniform and the magnification derived is very large, x200, 
independent of the source redshift. This model is used to 
identify 9 additional sets of fainter multiply lensed images, 
which are incorporated to improve the fit, so that the final 
model comprising 6 free parameters, is constrained by a to- 
tal of 33 lensed images. These fainter images are deflected 
by larger angles than the bright lensed spiral galaxy plac- 
ing their images generally further from the cluster center. 

From our model, the circular equivalent Einstein ra- 
dius contained within the critical area is ~ 27" for faint 
sources at z ~ 2, corresponding to the outer tangential 
critical curve (drawn in Figure 2), which in physical units 
is 170 ± 20kpc at the distance of the cluster. This very 
large radius adds to the already uncomfortable discrep- 
ancy between the large Einstein radii observed for mas- 
sive clusters and the predictions based on the standard 
LCDM cosmology (Broadhurst & Barkana, 2008, Sadeh & 
Rephaeli, 2008, Puchwein & Hilbert, 2009) for which such 
large Einstein radii can only be contemplated with mass 
distributions which are prolate and aligned along the line 
of sight (Corless k King, 2007, Oguri & Blandford, 2009). 
Instead, here the cluster is evidently elongated across the 
line of sight, traced by the bright galaxies, extended X-ray 
emission (see Figure 1 in Ebeling et al., 2007) and the crit- 
ical curves (see Figure 2). The central mass distribution 
is evidently unconcentrated, and presumably related to its 
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currently unrelaxed state. 

The velocity field of the lensed spiral can be readily ob- 
tained because of the large solid angle covered by each 
lensed image, and if the inclination can be reliably de- 
rived from such data then an independent estimation of 
the magnification can be derived via the TuUy-fisher re- 
lation. Spectroscopy of the internal velocity dispersion of 
the starlight of the cD galaxy will also help in understand- 
ing the cD mass profile, for which we have shown, rather 
sm-prisingly, Hc;cms to be dominated by stars out to the 
observable limit of 30 kpc. 

This cluster is unique in having near uniform density in 
projection, at the critical level, thereby maximising grav- 
itational lens magnification. We calculate that the total 
area of sky exceeding a magnification, fi > 10, is ~ 2.8 
square arcminutes, corresponding to the current high red- 
shift limit of 2; ~ 8, which is over twice the equivalent 
area calculated for other massive clusters such as A1689 



(Broadhurst et al., 2005) and C10024 (Zitrin ct al., 2009a). 
This extreme magnification together with the lack of im- 
age distortion makes MACS Jl 149.5-1-2223 the most pow- 
erful known lens for accessing faint galaxies in the early 
universe. 

ACKNOWLEDGMENTS 

We are grateful for the hospitality of ASIAA, where 
part of this work was accomplished. We thank our anony- 
mous referee for useful comments, and acknowledge use- 
ful discussions with Craig Sarazin, Holland Ford and Wei 
Zheng. ACS was developed imdcr NASA contract NAS 
5-32865. This research is based on observations pro- 
vided in the Hubble Legacy Archive which is a collab- 
oration between the Space Telescope Science Institute 
(STScI/NASA), the Space Telescope European Coordinat- 
ing Facility (ST-ECF/ESA) and the Canadian Astronomy 
Data Centre (CADC/NRC/CSA). 



REFERENCES 



Barkana, R., 1998, ApJ, 502, 531 

Bartclmann, M., Huss, A., Colberg, J.M., Jenkins, A., Pearce, F.R., 

1998, A&A, 330, 1 
Biviano, A., Salucci, P., 2006, A&A, 452, 75 
Bouwens, R.J., et al., 2004, ApJ, 616, 79 
Bradley, L.D., et al., 2008, ApJ, 678, 647 

Broadhurst, T.J., Taylor, A.N., Peacock, J.A., 1995, ApJ, 438, 49 
Broadhurst, T.. ct al., 2005, ApJ, 621, 53 
Broadhurst, T.J. & Barkana R., 2008, MNRAS, 390, 1647 
Broadhurst, T, Umctsu, K, Mcdczinski, E., Oguri,M., Repha«li, Y., 

2008, ApJ 685, L9 
Corless, V.L., King, L.J., 2007, MNRAS, 380, 149 
Dalai, N., Holder, G., Hennawi, J.F., 2004, ApJ, 609, 50 
Donnarumma, A., Ettori, S., Meneghetti, M., Moscardini, L., 2009, 

arXiv0902.4051 
Dressier, A., 1979, ApJ, 231, 659 

Ebbels, T.M.D., Le Borgnc, J.-F., Pello, R., Ellis, R.S., Kneib, J.-P., 

Small, I., Sanahuja, B., 1996, MNRAS, 281, 75 
Ebeling, H., et al., 2007, ApJ, 661, 33 

Franx, M., lUingworth, G.D., Kelson, D.D., van Dokkum, P.O. & 

Tran, K., 1997, ApJ, 486, 75 
Frye, B. & Broadhurst, T., 1998, ApJ, 499, 115 
Frye, B., Broadhurst, T., Bem'tez, N., 2002, ApJ, 568, 558 
Gavazzi, R., Fort, B., MeUier, Y., Pello, R., Dantel-Fort, M., 2003, 

A&A, 403, 11 



Gcbhardt, K., Thomas, J., 2009, arXiv0906.1492 
Halkola, A., et al., 2008, A&A, 481, 65 

Hennawi, J.F., Dalai, N., Bode, P., Ostrikcr, J. P., 2007, ApJ, 654, 
714 

Horesh, A., Ofek, E.O., Maoz, D., Bartelmann, M., Meneghetti, M., 

Rix, H.-W., 2005, ApJ, 633, 768 
Kneib, J.-P., Ellis, R.S., Santos, M.R., Richard, J., 2004, ApJ, 607, 

697 

Limousin, M., ct al., 2008, A&A, 489, 23 

Oguri, M., ct al., 2009, arXiv0901.4372 

Oguri, M., Blandford, R.D., 2009, MNRAS, 392, 930 

Puchwein, E. & Hilbert, S., 2009, arXiv0904.0253 

Sadeh, S. & Rephaeh, Y., 2008, MNRAS, 388, 1759 

Soucail, G., Fort, B., Mellier, Y., Picat, J. P., 1987, A&A, 172, 14 

Springel, V., et al., 2005, Nature, 435, 629 

Umetsu, K., et al., 2009, ApJ 694, 1643 

von der Linden, A., Best, P.N., Kauffmann, G., White, S.D.M., 2007, 

MNRAS, 379, 867 
Wambsganss, J., Bode, P., Ostrikcr, J. P., 2004, ApJ, 606L, 93 
Zheng, W., et al., 2009, ApJ, 697, 1907 
Zitrin, A., et al., 2009a, MNRAS, 396, 1985 

Zitrin, A., Broadhurst, T., Rephaeli, Y., Sadeh, S., 2009b, 
arXiv0907.4232 



